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Abstract 

We obtain new stronger bounds by orders of magnitude, on the ultimate 
temperature of the universe by exploiting the copious production of gravitinos 
at finite temperature. 
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In this paper we will derive new bounds on the temperature of the universe 
by considering the cosmological implication of gravitinos with masses varying 
from 10~ 4 eV to ITeV. The difference with previous results [1-6, 10] is due to 
the large interactions of the spin 1/2 component of the gravitinos with the hot 
primordial plasma. 

There are two major assumptions in this paper: 

1. We will assume that the theory describing ultramicroscopic physics leaves 
over at lower energies an N = 1 supergravity theory [7]. 

2. We will also assume that the typical splitting in supermultiplets after 

— 1/2 

supersymmetry breaking is of O (G F ). This represents the prejudice 
that supersymmetery "solves" the old fine tuning puzzle. 

The low energy spectrum then contains the usual particles and their super- 
partners as well as a gravitino. We will begin by discussing briefly the different 
masses for superpartners of usual particles. The origin for the mass of the 
gravitino is the superhiggs mechanism [8] which is of gravitational strength. In 
contrast, the masses of the other superpartners depend on how supersymmetry 
breaking is fed down to the low energy sector and may not involve to leading 
order the Planck mass, M p . This is why the gravitino mass 7713/2 can be quite 

— 1/2 

different than the typical superpartner mass of 0(Gp'') if the phy sics that 
mediate supersymmetry breaking to the usual sector, is not gravitational. Let 
me briefly review how masses of scalar partners and gauginos arise when grav- 
itational feeddown is negligible. The masses of gauginos come from dimension 
five operators of the form 



where 5 is a chiral multiplet whose F component breaks supersymmetry. M is 
the scale of the physics that mediates the breaking to the low energy sector. As 
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(1) 



an example [9] one could imagine heavy (mass M) superfields charged under 
the standard model gauge group SUc(3) x SUl(2) x Uy(l) that couple directly 
to the supersymmetry breaking sector. The low energy sector then feels the 
breaking of supersymmetry through usual gauge interaction with the heavy 
sector. In this case, M can be substantially lower than the Planck scale. Indeed, 
after supersymmetry breaking the dimension 5 operator leads to a gaugino 
major ana mass term 

where a are gauge indices. So that 

The scalar quarks, leptons and higgs masses appear through dimension 6 oper- 
ators of the form 

^2 / rf4 e S*S T*T, (4) 

where the T superfield stands for the quarks, leptons and Higgs superfields. 
After supersymmetry breaking, a (mass) 2 term is generated of the form 

\Fs\ 2 
M 2 

As stated before we will assume that 

\Fs\ 



A T At ■ (5) 



0(G F 1/2 ) . (6) 



M 

It is therefore possible to have F$ as small as a TeV and M of order 10 TeV. 
This would then lead to a tremendously small mass m 3 / 2 for the gravitino, 
m 3/2 ~ 0(lCT 4 eV). Indeed the gravitino mass due to the mixing with the 
would be goldstino is 

m 3/2 (7) 
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Then by using expression (6), m 3 / 2 can be reexpressed as 

1713/2 ~ ~M~ V ® 
Now in order for no fine tuning of scalar Higgs masses, m 3 / 2 should not be heav- 
ier than a TeV. In order to discuss the cosmological consequences of gravitinos 
we need to consider their interactions with other fields. As is well known [10], 
the lifetime r 3 / 2 for a gravitino of mass m 3 / 2 larger than the typical splitting 
AM in supermultiplets, is 

r 3/2 ~ M 2 p /ml /2 . (9) 

The lifetime r 3 / 2 then ranges from days to 10 4 years for m 3 / 2 in the range 
lOGeV <m 3 / 2 < ITeV. This presumes the lightest superpartner to be not 
heavier than lOGeV. As was pointed out in previous studies [1-6, 10], this 
could be potentially fatal for nucleosynthesis if gravitinos had earlier been in 
thermal equilibrium. Gravitinos lighter than the lightest superpartner could in 
principle decay to a neutrino and photon through higher dimensional operators 
that violate R parity and lepton number. An example is the dimension 6 
operator 

-L^*a»L a (F, u ) ab H b u (10) 

where ^ is the spin 3/2 component of the gravitino, L a is a lepton doublet, 
H\ is the Higgs doublet that gives mass to the up quarks and F^ is the SU{2) 
field strength. This operator leads to a lifetime r 3 / 2 for the gravitino 

r 3/2 r 5— (11) 

m 3/2 

which is much larger than the age of the universe for m 3 / 2 of O (G F 1/2 ^J . There- 
fore, such gravitinos are stable and the danger [1-6, 10], in this case, is that their 
present energy density, p 3 / 2 (today), could be larger than the critical energy 
density p c . 



The crucial information needed to discuss the cosmologial implications of 
gravitinos, is their abundance. In order to estimate this abundance, we will 
need to calculate their interaction rate with the early universe plasma. This 
is precisely where our analysis differs with earlier estimates [2-6]. When the 
energy E and momentum \k\ of the gravitino is larger than the typical splitting 
AM in a multiplet, but smaller than >J~Fs, the dominant interaction of the 
gravitino is through its longitudinal component, the "goldstino" Vl/g. In this 
range of energies one can derive the interaction of by using a non linear 
realization of supersymmetry [11]. Indeed, under a rigid supertransformation, 
the goldstino transforms as: 

5* a s = s a (F s ) + ... (12) 

so that it must decouple at zero momentum, as do goldstone bosons. The 
leading terms in the interaction Lagrangians are then: 



Lint ~ ps 



+{D p A)*d^sv P o^s ~ | (d^s) °» \*{A*T a A) 



dW 



+ h.c. (13) 



1 dF* 

where A a is a gaugino field, F* v is the field strength and A and \I> are respectively 
scalar and fermion components of chiral multiplets. One can then calculate for 
example in this regime the cross section a for gluon+gluon — ► gravitino+gluino, 
see fig. 1. As Fayet [12] pointed out, this cross section a is given by 

(AM) 2 , , 

a ~ a s ^yr~ ■ ( 14 ) 

Note that in the regime of energy AM < E < y/F, one could naively expect 
that due to the derivative coupling of the goldstino, the cross section depends 
quadratically on the energy: 

E 2 

a~a s —2. (15) 



This however does not occur, instead a saturates to the behavior found in 
expression (14). The reason for this, is that in this regime, the splitting AM 
becomes negligible compared to the energy scale so that supersymmetry is 
restored and the goldstino hence decouples. This implies that the leading energy 
dependence of the separate diagrams in fig. 1 cancels, leaving over (14). 

In the cosmological context however, we are interested in the cross section 
at high temperatures T in the range AM < T < y/F. Unlike the high energy, 
zero temperature case, the cross section a in this thermal environment does not 
saturate but instead grows as 



The reason for the growth in the cross section is that at finite temperature 
the bosonic states with momenta \k\ <T have large occupations, whereas the 
fermionic states are at most occupied by one quantum. At the leading order, 
depicted in fig. 1, because of energy-momentum conservation the cancellation 
still occurs [16]. However, the next order in as allows for intermediate bosonic 
states occupied by more than one quantum. Therefore, even though the tem- 
perature T is larger than the splitting AM, the leading temperature behavior 
at this next to leading order in as, will not cancel among the various diagrams. 
This is because the different statistics of bosons and fermions "break" super- 
symmetry in a thermal environment. The rate V for such interactions in this 
temperature regime by dimensional arguments is 




(16) 




rp5 



(17) 



The gravitino can therefore be in equilibrium at temperature 




(18) 



5 



If the temperature of the universe did reach such value, then the number density 
of gravitino n 3 / 2 (T) ~ T 3 . The energy density p^^iT) would then satisfy 

P3/2(T)~g*(T)m 3/2 T 3 (19) 

where g* is a dilution factor due to previous annihilations or decays, which 
implies for stable gravitinos p 3 / 2 (today) ~ 10 _2 to 3 / 2 (3 10~ 4 eV) 3 . We then 
recover the earlier literature result that gravitinos with masses m 3 / 2 < IkeV 
satisfy the constraint that their present energy density is bounded by p c . For 
masses m 3 / 2 in the range IkeV < m 3 / 2 < lOGeV, we find that these should 
never have been in equilibrium. Therefore, a simple upper bound on the maxi- 
mum Temperature T max is the Temperature T eq at which gravitino equilibrates, 

Te, ~ j~ s ™3/2 M P /3 ■ ( 2 0) 

There are presumably large numerical factors in (17), typical of thermal cal- 
culations involving large numbers of components, like in the Debye mass for 
gluons [17]. This is why the factor of ^ does not numerically affect much the 
bounds on T max . 

One can get a more serious bound by strictly implementing p 3 / 2 (today) 
< p c . This is done by integrating the relevant Boltzman equations. The result 
for T max as a function of m 3 / 2 is depicted on the diagram in fig. 2. Note that for 
m 3/2 ~ 10 GeV the maximum temperature T max ~ O (10 5 GeV). 

For unstable gravitinos with masses in the range 
lOGeV < M 3 / 2 ^ 1 TeV the bound on T max is obtained by requiring that nu- 
cleosysthesis not be affected. This entails that during this period, the gravitino 
doesn't dominate the energy density of the universe. Also, their decays after nu- 
cleosynthesis should not generate more entropy than was already present, such 
as not to dilute the abundance of light nuclei. The latter bound is stronger 
in constraining the abundance of gravitinos. One should also consider the de- 
struction of deuterium, and other light nuclei by decay products of gravitinos 
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[2, 3, 5, 6]. This was shown in the past, to lead to the most stringent bound 
on the ultimate temperature of the universe. 

The ratio of entropy before decay Si, to after decay Sf can be estimated to 

be 

S£ ^ K/2(73/2) ^3/2] 3/4 ^ 
Si "3/2 (^3/2) 

By equating the inverse lifetime of gravitinos with the expansion rate of the 

universe dominated by these gravitinos one obtains 

5 

"^3/2 

"3/2 (T-3/2) ~ -jp- • (22) 



This leads to 

1/2 



Si \m 3/2 J 

Requiring that the decaying gravitinos not modify in any significant way the 
already present entropy, leads to a bound on T max as a function of the gravitino 
mass m.3/2 as depicted on the diagram of fig. 3. For a range of gravitino mass, 
lOGeV ^m 3 / 2 < 1 TeV, the ultimate temperature T max range is then 

10 5 GeV<T ma:E <10 7 GeV. (24) 

If one considers instead, the bound from the destruction of light nuclei, following 
previous authors [2, 3, 5, 6] the result is even more dramatic in that it predicts 
a maximal temperature T max in the range 

10 4 GeV<T ma;E <10 5 GeV, (25) 

for the same range of gravitino mass, as can also be seen in fig. 3. These results 
are obtained by rescaling the bounds appearing in earlier papers [3,6], with the 
new temperature dependence of the gravitino abundance. 

In conclusion the possible existence of gravitinos with masses 
m 3/2 > 1 keV require low temperature baryogenesis. Such late generation 



of a baryon asymmetry was discussed earlier in a beautiful paper by Affleck 
and Dine [13]. Also, more recently, a flurry of activity around weak interaction 
baryogenesis [14] provides another way to reconcile observation with the possi- 
ble existence of gravitinos with masses of order or smaller than a TeV. Note in 
addition, that if the universe had experienced a period of inflation with cosmo- 
logical constant A, the reheating temperature Tr should then be smaller than 
Tmax- If h appears that the microwave anisotropy has a tensor component, 
then the cosmological constant A is of O[(10 17 GeV) 4 ] [15]. This would then 
imply an extremely inefficient reheating with 

Tr T max ,„„v 



A 1 / 4 10 17 GeV 

Tr 

Therefore the ratio should not exceed an extremely small number in 

the vicinity of 10~ 12 , 10 3 . This in turn requires an extremely weak coupling 
between the sector responsible for inflation and the standard model degrees of 
freedom. This estimate is conservative since the production of gravitinos by 
the inflationary sector during the exit of inflation was neglected. 

To avoid any bound on the temperature of the universe, gravitinos should be 
lighter than 1 keV. This in turn implies low energy supersymmetry breaking at 
a scale F$ comparable or smaller than (10 6 GeV) 2 . In this case, supersymmetry 
has to be mediated to the low energy sector by forces other than gravity. 
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Figure Captions 

Fig. 1: Feyman diagram is involved in the calculation for the cross section of 
gluon + gluon — ► gluino + "goldstino". 

Fig. 2: The diagram for the ultimate temperature T max as a function of m 3 / 2 , in 
the case of stable gravitinos. 

Fig. 3: The dia gram for the ultimate temperature T max as a function of m 3 / 2 in 
the case of unstable gravitinos. 
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This figure "figl-l.png" is available in "png" format from: 
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Figure 1 
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This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-th/9404044v4 
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This figure "figl-3.png" is available in "png" format from: 
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Figure 3 
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